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Harvesting Excitons through Plasmonic Strong Coupling 
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Exciton harvesting is demonstrated in an ensemble of quantum emitters coupled to localized surface 
plasmons. When the interaction between emitters and the dipole mode of a metallic nanosphere reaches the 
strong coupling regime, the exciton conductance is greatly increased. The spatial map of the conductance 
matches the plasmon field intensity profile, which indicates that transport properties can be tuned by 
adequately tailoring the field of the plasmonic resonance. Under strong coupling, we find that pure 
dephasing can have detrimental or beneficial effects on the conductance, depending on the effective number 
of participating emitters. Finally, we show that the exciton transport in the strong coupling regime occurs 
on an ultrafast timescale given by the inverse Rabi splitting (~ 10 fs), orders of magnitude faster than 
transport through direct hopping between the emitters. 

PACS numbers: 71.35.-y, 05.60.Gg, 73.20.Mf, 81.05.Fb 


Bound electron-hole pairs in semiconductors and molec¬ 
ular solids, known as excitons, play a key role in many 
basic processes such as Forster resonance energy transfer 
or energy conversion in light-harvesting complexes [1-3]. 
Various optoelectronic devices are also based on exciton 
dynamics, including organic solar cells [4], light-emitting 
diodes [5], and excitonic transistors [ 6 ]. Thus, a large re¬ 
search effort is directed towards controlling exciton trans¬ 
port properties. As excitons usually suffer from relatively 
large propagation losses associated with decoherence and 
recombination, increasing their propagation length is an 
important goal. Moreover, the transport of these quasiparti¬ 
cles into a specific spatial region in a controlled way could 
significantly increase the efficiency of devices such as solar 
cells, where the presence of excitons close to the charge 
separation region is essential for photocurrent generation. 

Recent works [7, 8 ] have shown that exciton conductance 
in an ensemble of organic molecules can be boosted by 
several orders of magnitude when the ensemble is coupled 
to a cavity mode and the system enters the strong coupling 
(SC) regime. This is a promising result, which motivates 
the exploration of this phenomenon beyond cavity setups. 
Among the fields in which the SC regime is very attractive, 
plasmonics stands out due to the recent works studying SC 
between surface plasmons and quantum emitters (QEs) [3, 9- 
12, 14-20]. Surface plasmons arise as ideal platforms for 
SC applications, due to their small mode volume and the 
tunability of their electric field profile via nanostructure 
design. These characteristics could allow for a deterministic 
control of exciton harvesting. 

In this Letter we demonstrate efficient harvesting of ex¬ 
citons in a collection of QEs strongly coupled to localized 
surface plasmons (LSPs). As a proof of principle, we first 
study a system of QEs interacting with the dipolar modes 
supported by a metal nanosphere (NS). We show how within 
the SC regime the spatial map of the exciton conductance 
mirrors the field intensity profile. Employing a structure 
composed of three aligned nanospheres, we show that the ef¬ 


ficiency of exciton harvesting can be significanly increased 
by tuning the electric field profile of the LSP mode. In 
addition, we demonstrate that the role of dephasing in the 
exciton conductance strongly depends on the effective num¬ 
ber of emitters involved in the formation of the polariton 
modes. Finally, we demonstrate that the speed of exciton 
transport in the SC regime is orders of magnitude faster than 
in the weak coupling regime. 

The first system we consider consists of a silver NS of 
radius R, surrounded by a layer of N QEs, as shown in 
Fig. la. For simplicity, the QEs are regularly distributed 
over a spherical layer, which we place at a distance h — 
1 nm away from the surface of the NS. Note that at this short 
distance, higher multipole modes of the NS are dominant 
and lead to quenching losses for a single QE [21]. However, 
recent works have shown that for N emitters, collective SC 
with the dipole resonance of the NS does indeed arise, and 
higher multipoles merely add an effective detuning to the 
hybrid mode [18]. Hence, for the silver nanosphere we only 
consider the three dipolar LSP modes (x, z), which are 

characterized by their frequency oUpi and decay rate tz, and 
electric field profile Eo,{r) {a G x, y, z). These parameters 
can be extracted from the NS properties, as described in the 
Supplemental Material [22]. The QEs are modelled as point 
dipoles oriented along the radial direction, with transition 
frequency ouq, dipole moment jl, and total decoherence 
rate 7 = 7 ^ + 7 d, where 7 ^ accounts for pure dephasing, 
while the decay rate 7d = 7r + 7nr contains radiative and 
nonradiative contributions. 

The Hamiltonian associated with the emitters-NS system 
within the rotating wave approximation can be written as 

H='^uioc]cj+ uipialaa+ 

j a=x,y,z 

X + '}2^9jac]aa + H.C.). 

iAj 

Here, the operators Cj and annihilate an excitation in 
emitter j (= 1,..., A^) and the LSP mode a, respectively. 
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The dipole-dipole interaction between QEs is given by Vij, 

and the QE-LSP coupling is Qjo, — —jl • Ea{rj). As 
we will show later, in the SC regime only a single ESP 
mode contributes and dipole-dipole interactions can be ne¬ 
glected. Under these approximations and for zero detuning 
(ojq = (jjpi), the + 1 singly excited eigenstates of H 
are formed by: i) two polaritons |zb) = ^(a^|0) zb \B)), 

where \B) = ^ the collective molecular 

bright state, with the Rabi splitting (12^ = 4 P)’ 

and ii) the so-called dark states, — 1 combinations of 
molecular excitations orthogonal to \B) which have no ESP 
component. The eigenfrequencies of the two polaritons are 

CUq zb ^r/2. 

In order to study exciton transport through the ensem¬ 
ble of emitters, we first determine the steady state of 
the system when one of the QEs (emitter A from now 
on) is incoherently pumped. Notice that this is the only 
driving term and no additional external illumination is 
present. The system is described by its density operator 
p, whose dynamics is governed by the master equation 
P = -i[H,p] + InCO- 

herent processes (losses and pump) are described by Eind- 
blad terms, Ci,[p] — bpb^ — {b^b^ p}/2. In the first part 
of this work we treat pure dephasing as a decay chan¬ 
nel for simplicity. The pumping rate 7^ is chosen small 
enough to stay in the linear regime. The steady state density 
matrix pss is obtained numerically using the open-source 
QuTiP package [23]. Einally, the exciton conductance mea¬ 
suring the efficiency of exciton transfer from emitter A 
to emitter j is calculated as = J(^A^j)l'yp, where 
J{A^j) = 7 Tr {HCcj [Pss]) is the energy loss rate of emit¬ 
ter j [7]. 

The exciton conductance is shown in Eig. lb for QEs at 
three representative positions, B, C, and D, as depicted in 
Eig. la. With organic molecule applications in mind, the 
parameters of the W = 100 QEs are chosen to correspond 
to TDBC J-aggregates at room temperature [24-26]: ouq — 
2.11 eV, p — 0.749 e-nm, 7 ^ = 26.3 meV, 7 ^ = 1.32 • 
10“^ eV, 7nr = 1.10 meV. The nanosphere (radius R — 
10 nm) is embedded in a dielectric host of permittivity = 
6 . 8 , so that the ESPs are in resonance with the QEs. Einally, 
the plasmon losses are given by = 0.1 eV. In order to 
study different Rabi frequencies we first artificially 
tune the field strength of the ESPs, instead of varying the 
number of QEs as could be done experimentally. Our results 
show how the onset of SC clearly differentiates two different 
regimes for exciton transport. In the weak coupling case 
(small Qi^), the dipole-dipole interaction between the QEs 
governs the dynamics and the transport is rather inefficient 
over large distances. The exciton conductance to emitter 
D is thus smaller than that to B or C. In the SC regime, 
however, this situation changes drastically, as ESP-mediated 
interaction becomes the primary transport channel [7]. Due 
to the dipolar field profiles, emitter A only couples to the 
z—dipole mode of the NS, and the x and y dipole play a 
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FIG. 1. (color online) a) Illustration of the single-NS system. The 
colored background shows the electric field intensity associated 
with the z-oriented dipole mode, b) Exciton conductance versus 
Rabi splitting for = 100 emitters. The gray dashed line indi¬ 
cates the onset of SC, Qr > I 7 — /^|/ 2 . c) Angular dependence 
of the conductance in the SC regime (light green) and the same 
quantity when neglecting dipole-dipole interaction between QEs 
(dark green). For comparison the three-NS case is also depicted 
(black line). 


negligible part. As the couplings gja are proportional to 
the electric field of mode a, excitons are transferred more 
efficiently to regions of high field intensity. This is the cause 
of the boost in the conductance of emitter D displayed in 
Fig. lb. 

These results indicate that the strong coupling conduc¬ 
tance is position-dependent, mimicking the field intensity 
profile. This fact is confirmed in Fig. Ic. Here, we plot the 
conductance in the SC regime for every QE versus its po¬ 
lar angular coordinate 6. The results without dipole-dipole 
coupling (Yij — 0) agree very well with the full calculation. 
The significant dip of (jg around 9 — 90° confirms that only 
the z-oriented ESP mode plays a relevant role. Therefore, 
we can safely neglect both the dipole-dipole interactions 
and the x- and ^-oriented ESP modes within the SC regime. 
Under these approximations, an analytical solution for the 
master equation can be obtained [22]. For zero detuning 
((jjpi — (jJq), the exciton conductance between emitters A 









and j has the simple form 
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0) ^ 16|g^ngjfcJo (7 + r) 

r(fi| + 7K)(fi| + 27 r)’ 

where we have defined the rate F = 7 + and Qj = for 

simplicity. Up to now, we have considered the case in which 
both the pumping and collection involve a single QE, and 
Q.R is artificially modified by tuning the field strength of the 
ESP while keeping N constant. In a realistic experiment, 
several emitters near location A would be pumped, and 
excitons collected from a region around D. Furthermore, 
oc y/N would be varied by changing the number of 
emitters N. Equation (2) is easily generalized to this case, 
giving 


VAriD(^o{l + ^)^% 

r(f2U7«^)(fiU27r)’ 

where rjx = measures the fraction of 

the Rabi frequency due to emitters involved in the pumping 
{tia) and collection processes {rjo), respectively. Both rjA 
and rjD are independent of N for uniform distributions 
of emitters. For small N (i.e., weak coupling), Eq. (3) 
shows that (jg grows as N‘^, while when the SC regime is 
entered for large N 7 r), it saturates to a constant 

value, r]AT]D(^o{^ + "f/^)> This equation thus predicts that 
the exciton transport efficiency from a pumping site to a 
collection spot can be increased by tailoring the mode to 
have maximal field strength at (only) these two locations in 
order to have large tja and 7]^. 

We demonstrate this by adding two additional identical 
silver nanospheres to the existing structure, as shown in 
Fig. 2a. The background of the panel displays the field 
intensity map of the lowest energy mode, which in this case 
is not degenerate as the rotational symmetry is broken. The 
spheres are separated by a 2 nm gap. In order to facilitate 
the comparison with the single-NS case and focus on the 
effect of the different mode profile, the ESP frequency and 
losses as well as the QE properties and locations are kept 
unchanged. Furthermore, we again plot the single-emitter to 
single-emitter conductance. The exciton conductance as a 
function of the Rabi splitting is displayed in Fig. 2b for N = 
100. As the figure shows, the onset of SC again produces 
a substantial increase in the pole-to-pole conductance 
significantly larger than in the single-NS case. The most 
striking feature of this system is that transport to emitter D 
is now much more efficient than to the nearest neighbor B 
of emitter A. The conductance in the SC regime is much 
more concentrated around the poles than in the single-NS 
(see Fig. Ic), and conductance to point D is increased by a 
factor of 50. Excitons are thus shown to be very efficiently 
harvested at the hot spots of the ESP mode. This is an 
interesting result towards potential applications, due to the 
high tunability provided by the wide variety of plasmonic 
nanostructures that are available nowadays. 



^i?(eV) 


FIG. 2. (color online) a) Illustration of the three-NS system. The 
colored background shows the electric field intensity of the lowest 
energy mode, b) Exciton conductance versus Rabi splitting for 
A = 100 emitters. The curves correspond to transport from 
emitter A to emitters B, C, and D as depicted in panel (a). As in 
Fig. lb, the gray dashed line indicates the onset of SC. 


Up to now, our treatment of relaxation processes has been 
relatively crude, lumping together dissipation and dephas¬ 
ing. Recent works show that dephasing mechanisms can 
be relevant for exciton transport in organic compounds [2]. 
Therefore, we now describe pure dephasing in more detail 
by performing the substitution 7 ^^ ^ Id^cj + 
in the master equation. In order to avoid unrealistic conclu¬ 
sions, we have also checked that the complementary Bloch- 
Redfield-Wangsness formalism [27] reproduces the general 
behavior obtained with the standard Eindblad method. 

The pole-to-pole exciton conductance in the SC regime 
is shown in Fig. 3 as a function of the dephasing rate, for 
both the single-NS and the three-NS cases. Except for the 
value of 70 , the parameters of the three nanospheres and the 
QEs are the same as in the previous calculations. Surpris¬ 
ingly, the dependence of the conductance with dephasing 
is remarkably different for the two considered nanostruc¬ 
tures. As dephasing is increased, the conductance decreases 
monotonically for the three-NS structure, but counterintu¬ 
itively increases for the single NS. As we show next, this 
difference in behavior is related to the effective number of 
QEs that enter strong coupling with the ESR This number 
is quite small in the three-NS case, where only QEs near the 
hot spots participate in SC. Our hypothesis is confirmed by 
considering a system of N non-interacting QEs uniformly 
coupled to a cavity mode (pj — g). In this simple situation, 
an analytical formula for the exciton conductance can be 
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FIG. 3. Exciton conductance from point A to point D in the 
SC regime = 1 eV), as a function of the dephasing rate. 
The single-NS case (orange line) is compared with the three-NS 
structure (red line). The inset displays the conductance in the SC 
regime for a system of QEs uniformly coupled to the field, in the 
cases N = b and N = 60. 

derived [7], which in the SC limit is given by 

^ ^ ^o7d(7d + + 27d + 270 ) 4 

(270 + 7dAr) + 7diV(K + 7d + 7^)) ■ 

This simple expression (shown in the inset of Fig. 3) is 
able to reproduce the observed features. Specifically, the 
ratio (Je{l(j) C)c)/(Je (70 ^ 0 ) is approximately equal to 
Since typical plasmonic structures fulfill Ad/^ 
1 , this expression predicts that with increasing dephasing, 
exciton conductance decreases for small N (N < \//Ad)^ 
but increases for sufficiently large N (N > \/i^/Ad)- This 
dependence with the number of emitters suggests that the 
dark states play a key role in this process. Since dephasing 
creates an incoherent coupling between the dark modes 
and the polaritons, a fraction of the population in the dark 
modes can be transferred to the polaritons. For large N, the 
dark states are highly populated as the overlap between the 
initial state (one excitation at emitter A) and the polaritons is 
extremely small. As a consequence, the excitation transfer 
to the polaritons can compensate for the detrimental effect 
of dephasing on this state, making the conductance increase. 

While we have so far focused on the conductance as 
obtained in the steady state under pumping, the temporal dy¬ 
namics of the system provides important additional insight. 
We thus investigate the population dynamics in the single 
NS case, for an initial excitation of emitter A. In the weak 
coupling regime, the dipole-dipole interaction dominates 
and slowly transfers population to emitter D (as shown in 
Fig. 4a). The plasmon modes do not significantly participate 
in the dynamics, so that all populations decay with the life¬ 
time of the bare QEs (r ^ 600 fs) for large times. As the 
exciton transport is even slower, the increase of population 



FIG. 4. (color online). Time dynamics of the single NS system 
when QE A is initially excited (initial state c^|0)) and no pump is 
applied. The populations of both the main ESP and the relevant 
emitters, A and D, are displayed as a function of time, a) Weak 
coupling case, = 10“^ eV. b) Strong coupling case, Qji = 
0.3 eV. In both panels, the linear scale for short times allows for a 
better visualization of the results. 

in emitter D is cut off at around this time, with a maximum 
population of ^ 2 • 10“^. On the other hand, in the SC 
regime (Fig. 4b), the population is delivered to the emitter D 
much more efficiently through Rabi oscillations. These pro¬ 
ceed on a timescale determined by the inverse of the Rabi 
splitting (1/^R ^ 15 fs in Fig. 4b), giving extremely fast 
population transfer. Furthermore, the population of emitter 
D reaches significantly larger values than in the weak cou¬ 
pling regime, up to ^ 4 • 10“^. For large times, most of 
the population is trapped in the dark states, which can now 
also decay through the dephasing-induced coupling to the 
polaritons, giving an effective lifetime somewhat below the 
bare QEs. 

In conclusion, we have demonstrated the feasibility of 
harvesting excitons through strong coupling in systems com¬ 
posed of a plasmonic structure and an ensemble of organic 
molecules. We have shown that for emitters coupled to the 
localized surface plasmons of a single nanoparticle, the exci¬ 
ton conductance map mimics the electric field profile of the 
plasmon resonance. Taking advantage of this property, we 
have devised a more complex structure in which an exciton 
can be efficiently transferred between two subwavelength 
hot-spots of the plasmonic system. We have also shown 
how dephasing can be beneficial or detrimental depending 
on the number of emitters that are effectively coupled to the 
plasmon resonance. We have additionally demonstrated that 
exciton transport in the strong coupling regime proceeds or¬ 
ders of magnitude faster than under weak coupling. Finally, 
it is worth noting that our findings regarding harvesting of 
excitons mediated by strong coupling are general and appli- 
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cable to any quantum emitter, from atoms and quantum dots 
to organic molecules, and also to any confined electromag¬ 
netic mode with similar properties as the plasmonic ones 
used here. 
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Supplemental Material 


CALCULATION OF THE MODE PROPERTIES 


We calculate the classical field profiles and modal charac¬ 
teristics of the localized surface plasmon (LSP) numerically 
with the finite element method (using COMSOL Multi¬ 
physics) to solve Maxwell’s equations. The permittivity of 
the silver structures is given by a Drude-Lorentz formula: 

/X A 

e(w) = Coo - r^- —V “ — 02 I • r ’ 

where the parameters Co© = 3.91, = 8.833 eV, — 

0.0553 eV, A = 0.76, Vtp = 4.522 eV, and Vp = 8.12 eV 
are taken from Ref. [1]. For the dipole mode of the single 
NS, we obtain a frequency uj^i = 2.11 eV and a dissipation 
rate = 0.12 eV. In the three-NS structure, the parame¬ 
ters of the lowest energy mode are — 1.85 eV, and 
Kj' = 0.14 eV. Note that in the main text these values are 
modified artificially to agree with the single-NS case, in 
order to isolate the influence of the electric field profile. 


QUANTIZATION OF LSP FIELDS 

The classical fields of the system modes obtained from 
Maxwell’s equations have to be properly quantized in order 
to include them in a quantum Hamiltonian. In this section 
we briefly detail the method followed for this purpose. 


Single-sphere case 


The dipole mode of a single sphere is quantized by com¬ 
paring the classical and quantum values of the polarizability 
a. For a small metallic sphere of permittivity e(cj) and 
radius i?, the classical static polarizability is given by [? ] 


a cl = 


e(^) - Q 

e(w) + 2erf ’ 


( 2 ) 


being the permittivity of the surrounding dielectric. For 
simplicity, we introduce an approximate Drude permittivity 
(Eq. (1), with A = 0). We can define a resonance frequency 
ujr = uOp /+ 26^ and, after some algebra [3], arrive to 
the following expression: 


ad ^ -27TeoedR^ 


3ew 


cc;: 


+ 2ed cj(cj — uJr -\- i'yDl2) 


(3) 


which is valid in the vicinity of a narrow resonance, i.e. 
U ^ Ur ^ 7P)- 

On the other hand, the polarizability of a quantum two- 
level system with dipole moment fi is given by [4]: 

fR 2u() 


Uq — (u + 


2 ’ 


(4) 


where Uq and 7 are the two-level system frequency and 
linewidth, respectively. Close to a narrow resonance (u ^ 
^0 ^ 70)5 the above expression is approximated as 


^_CJo_ 

h u {u — Uq -\- i 7 o/ 2 ) 


(5) 


By direct comparison of Eqs. (3) and (5), we obtain the 
dipole moment of a nanosphere in a quantum model: 





Coo + ‘2.^d 


1/2 


( 6 ) 


The quantum electric field in this case corresponds to the 
classical field emitted by a dipole with dipole moment fipp 


Arbitrary geometry 


For more complex nanostructures, the eigenmodes do 
not possess a purely dipolar profile, and a new quantization 
procedure is necessary. We start with the classical electric 
and magnetic field profiles obtained from our simulations, 
Ed{r) and Hd{E). We generalize the usual quantization 
procedure for vacuum fields [4], by applying the correspon¬ 
dence principle to the electromagnetic energy: 

u(^Eci,Hci) (7) 

In the above expression, U stands for the total classical 
electromagnetic energy, and H for the Hamiltonian operator. 
The quantum field operators Eq^ Hq are related to their 
classical counterparts, 

P^CEcia + C*E:,a\ 

H, = CHcia + C*H:ia\ 


where a and are the photonic mode annihilation and 
creation operators, respectively, and C is a normalization 
constant we have to determine. 

Our systems are composed of a metallic nanostructure 
surrounded by a dielectric host, with permittivities e{u) 
and Cd respectively. Usually, the electromagnetic energy 
U is ill-defined in lossy media, and a macroscopic QED 
formalism is required. However, in the vicinity of a nar¬ 
row resonance u = u'^i, and provided that the losses are 
small (Im[e((x;)] <C Re[e(cc;)]), the following approxima¬ 
tion holds [5] 


U K — 
2 


+ 




L 

I 

7Vdiel 


dVRe 


r^metal 

dVE, 


d{ue{u)) 


cl 


du 


e; . E 


cl 


I 

JVt 


( 9 ) 

dVHii • iici. 


2 
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We apply the correspondence principle, Eq. (7), by intro¬ 
ducing the quantum fields (8) in the above equation, i.e. 
substituting 0*1, E^i, H*^, H^i) by {E]^,Eq, Hq). Af- 
ter manipulation, we arrive to the familiar Hamiltonian 


Tavis-Cummings expression, 

N 

H = 6^0 ^ ^ CdplCi^ Ci + (90^3 + H-C.) • (12) 

i=i j 


H = |C|V(2a^a + l). 


It is useful to work in the bright-dark basis, which is com- 
^^posed of a bright state \B), 


The comparison with the Hamiltonian of a harmonic oscil¬ 
lator, H = h(jj'pi{a^a + 1/2), gives the expression of the 
normalization constant. 





( 11 ) 


Finally, after C is determined, the calculation of the cou¬ 
plings g'j = —jl-C- Eci{fj) is straightforward. 

It is important to mention that, when calculating the en¬ 
ergy U from our simulations, we have to deal with the 
linear divergence of the second integral in Eq. (9). This is 
a fundamental problem regarding lossy cavities, in which 
normal modes are ill-defined. In our case, we can safely 
ignore the divergent contribution due to the low loss rate [6]. 
The validity of this approximation has been checked in the 
single nanosphere case, where this method reproduces the 
analytical results obtained above very accurately. 


ANALYTICAL FORMULA FOR THE CONDUCTANCE 


Here we give details of the calculation of the formula for 
the conductance (Eq. (2) in the main text). In the absence of 
dipole-dipole interaction, and considering only one of the 
dipole modes, the Hamiltonian of the system is given by the 


2 ^ 

l^) = ?rE5ic]|o), (13) 

“«i=i 

and a set of (W — 1) dark states. In principle, we are free 
to choose any orthonormal set of states, and we use this 
freedom to our advantage. As in our problem we pump the 
first molecule (state c\ |0)), we choose the dark states such 
that this pumping only excites one of them, which we name 
\D). It is straightforward to prove that 

10 ). (-) 

where we have defined G'^ = — g\ for simplicity. 

The remaining dark states {|i9/c); k = 1, ...W — 2} are 

arbitrary, provided that they fulfill the orthogonality condi¬ 
tions {Dk\B) = {Dk\D) = 0. Note that Eqs. (13) and (14) 
imply that {Dk\(jl\0) = {Dk\{a\B) + b\D)) = 0. 

Consequently, the states {D]^\ are not coupled to 
the pumped state and do not take part in the dynam¬ 
ics. The time evolution of the system is thus restricted 
to the 4-dimensional subspace spanned by the states 
{|0), |S), |79), |0)}, and the master equation reduces to 

a 16 X 16 linear system. 

We can thus calculate the steady-state density matrix 
analytically, from which it is straightforward to obtain the 
exciton conductance. With the notation employed in the 
main text, the general expression is given by 


| 2 | + 47 r^) + Up 0%r{T + 7) - 27(27^ + 5^7^ - + 47K^ + K0 

^e(J) \ 9 j \ \ 9 a \ (-f) 4 +47(452 + r 2 ))(f 22 jr 2 +7^(452 +r 2 )) 


where we have defined the detuning S = cjsp — cjq. The 
case of zero detuning 5^0 reduces the above expression 
to Eq. (2) in the main text. 
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